Multipartite Leggett-type Inequalities 
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We use two different approaches to derive multipartite Leggett-type inequalities, which are gen- 
eralizations of the two-qubit Leggett-type inequality obtained in [Nature Phys. 4, 681 (2008)]. The 
first approach is based on the assumption that the probability distributions should be non-negative. 
The second approach is based on a very simple algebraic equation and is, to some extent, easier 
than the first approach. Although these inequalities might not be the optimal ones in the sense that 
their quantum violations may not be the strongest, our results make the first step of generalizing 
Leggett-type inequality to multi-qubit systems and provide the possibility to experimentally test 
non-local realism in such systems. Moreover, the two approaches here may shed new light on the 
challenging problem of obtaining stronger multipartite Leggett-type inequalities. 
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The concept of physical realism suggests that the re- 
sults of observations are determined by the intrinsic prop- 
erties of a physical system and independent of measure- 
ments |l[ . This concept has taken root in classical physics 
and its significance goes far beyond science. Quantum 
mechanics (QM), however, challenges this concept in a 
very deep way. In 1964, Bell published his celebrated 
inequality based on Einstein-Podolsky-Rosen's notion of 
local realism (LR) [||. It was shown that any local re- 
alistic theory should obey Bell's inequality, while it can 
be violated easily in QM. Today this brilliant inequality 
has had various generations and many experiments 
have been carried out to test these Bell's inequalities [|| . 
Notwithstanding some loophole problems ||, these ex- 
periments overwhelmingly support QM and show viola- 
tions of Bell inequalities, and thus rendering local hidden 
variable models untenable. 

Then, should non-local realism hold water? To an- 
swer this fundamental question, Leggett in 2003 made 
a significant step by proposing an alternative non-local 
hidden variable (NLHV) model that was proved to be at 
variance with quantum predictions [10 . He introduced 



a new inequality based on this model and showed that 
this inequality can be violated by quantum correlations. 
Recently, a series of experiments have been carried out to 
test such model pd|-[T^|. These experiments again favor 
QM, casting doubt on the validity of non-local realism. 
Nevertheless, all these works only concern with two-qubit 
systems. Up to now, we still lack multipartite Leggett- 
type inequalities that can be used to experimentally test 
nonlocal realism in multi-qubit systems. 



In this paper, we generalize the Leggett-type inequality 
to multi-qubit systems by using two different approaches. 
The first approach follows a recent work by Branciard 
and his collaborators [^3|. While the second approach 
is based on a very simple equation and is a little bit 
easier. These inequalities are all satisfied by Leggett's 
NLHV model and violated by QM. Thus, our results have 
paved the way for experimental test of nonlocal realism 
in multi-qubit systems. However, a drawback of these 
inequalities is that the quantum violations do not scale 
with the system size as the violations of Bell inequalities 
do H 15 . In fact, the maximal quantum violations are 
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the same as for two-qubit case. 

For convenience and simplicity, we first focus at three- 
qubit system and the generalization to n-qubit case is ob- 
vious. Consider a common experimental scenario: three 
observers, denoted by A, B, and C (or Alice, Bob and 
Charlie), perform experiments with settings labeled by 
a, b and c, respectively. Their outcomes are denoted by 
a, f3 and 7 (a, /3, 7 = ±1). For the qubit case, a, b and c 
are vectors on the Poincare sphere and are independently 
and freely chosen by Alice, Bob and Charlie. Originally, 
Leggett's NLHV model is based on pairs of photons and 
three main assumptions [|l0| pT] : (i) Realism. All mea- 
surement outcomes are predetermined and independent 
of the measurements, (ii) Definite polarization. Physical 
states are statistical mixtures of subensembles. Within 
each subcnscmble, every photon in the pair has definite 
polarization, (iii) Malus's law. Within each subcnscm- 
ble, local marginals should obey Malus's law. Based on 
these assumptions, Leggett was able to derive an inequal- 
ity that was shown to be violated in QM. In our simpli- 
fied formulism, assumption (i) means that the outcomes 
of each observable is predetermined by some set of hid- 
den variables A, polarization parameters u, v and s (u, v 
and s are all unit vectors), and some set of other possible 



2 



non-local parameters x (here for simplicity, we choose 
these non-local parameters to be measurement settings 
in space- like separated regions O] ) . Mathematically, we 
have a = a(A, u, v, s, a, b, c), /3 = /3(A, u, v, s, a, b, c) 
and 7 = 7(A, u, v, s, a, b, c). Note that here the local- 
ity requirement in local hidden variable theory is explic- 
itly removed. This implies a big difference in deriving 
Bell inequalities and Leggett-type inequalities. For as- 
sumption (ii) and (hi), we denote the polarization dis- 
tribution function of subensembles and the probability 
distribution of A in each subensemble by D(u, v, s) and 
Pu,v.s(A), respectively. Then the expectation values of 
measurements for each subensemble are given by average 
of the measurement outcomes over the probability dis- 
tribution /9 u ,v.s(A). For example, we have ct/3(u,v,s) = 
/ dA/3 u ,v,s(A)a(A, a, b, c)/3(A, a, b, c) and a/3j(u,v,s) = 
J d\p UyV , s (\)a(\, a, b, c)/3(A, a, b, c)7(A, a, b, c). Fur- 
thermore, the assumption of Malus's law indicates 
a(u) = J dA/9 UiV ,s(A)a(A, u, v, s, a, b, c) = u • a, and sim- 
ilarly /3(v) = v • b, 7(s) = s • c. The assumption (ii) also 
implies that the physically measurable correlation func- 
tions are given by averaging the expectation values over 
the subensemble distribution D(u, v, s). For instance, the 
three body correlation function is given by: 



can be respectively expressed as: 



Q = (afa) 



J dudvdsD(u, v, s)aj3j(u, v, s) 



(1) 



After the introduction of the three major assumptions 
and the basic notions, now let us derive the three-qubit 
Leggett-type inequality using two approaches. The first 
approach follows the recent result of Ref. |13f . To this 
end, note that the conditional probability distribution 
P(a, /?, 7|a, b, c) can be expressed as: 



P(a, /3, 7|a, b, c) = / dKD(n)V K (a, /3, 7 |a, b, c), 



(2) 



where we use a single parameter n to denote (u, v, s) 
and P K (a, /3, 7|a, b, c) is the conditional probability 
of subensemble labeled by n. In the qubit case, 
P K (a, (3, 7|a, b, c) can be decomposed as: 

P K (a,/3, 7 |a,b,c) = i[I + c<(a,b,c)+/3£ B (a,b,c) 

+ 7 £ c K (a,b,c)+a/3£f(a,b,c) 
+ a7 £t C (a,b,c)+/37^ c (a,b,c) 
+ a/37C C (a,b,c)] (3) 

One advantage of this expression is that it enables one to 
clearly distinguish the marginals and the correlation co- 
efficient as discussed in Ref. |l3) . From the above Eq. (||) , 
it is obvious that £*(a, b,c), £ B (a, b, c), £^(a, b, c), 
£f(a,b, c), £f(a,b,c), £ BC (a,b,c), and £™ (a, b, c) 
have their physical meaning of average within subensem- 
ble k. For instance, Alice, Bob and Charlie's marginals 



£ B (a,b,c) = /3P«(a,/3,7|a,b,c) = j8(v), 
£°(a,b,c) = ^ 7P K (a,/3,7|a,b,c) =7(s). 

Similarly, the two-qubit correlation coefficients read: 
£t B (a,b,c) = a/3P«(a,/3,7|a,b,c) = 5ft 

q,/3,7 

£* C (a,b,c) = ^ a7P K (a,^,7|a,b,c) = 07, 

a,/3,7 

£ BC (a,b,c) = ^ /? 7 P K (a,/3,7|a,b,c) 

a,/3,7 

and the three-qubit correlation coefficient reads: 
£f c (a,b,c) - J2 «/^7P K («,A7|a,b,c) = 5£t~. 

a,/3,7 

As in Ref. |l3|, we also only concentrate on correlations 
satisfying the so called no-signaling condition, which im- 
plies the independence of marginals on other observer's 
inputs. For the marginals on Alice, Bob or Charlie's side, 
the condition is already certified by Malus's law O]. The 
no-signaling condition also indicates three more equa- 
tions: Cf (a,b,c) = £f (a,b), Cf (a,b,c) = £* c (a,c), 
and £ BC (a,b,c) = £ BC (b,c). The probability distri- 
butions P K (a, (3, 7|a, b, c) should be non-negative for all 
a, /3,7= ±1. This leads to eight inequalities: 
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0, (4a) 

0, (4b) 

0, (4c) 

0, (4d) 

0, (4c) 

0, (4f) 

0, (4g) 

0, (4h) 



From the inequalities (Ea-4h), it is easy to obtain: 

|£ A K (a)±£ BC (b,c)|<l±£r(a,b,c). (5) 

Now let us consider the case that Alice has two measure- 
ment settings a and a', while Bob and Charlie have only 
one measurement setting b and c, respectively. Then we 
can obtain four inequalities from inequality (H), two for 
(a, b, c) and two for (a',b,c). Combining these four in- 
equalities and using the triangle inequality, one arrives 
at: 



|£* BC (a,b,c)±£A BC (a',b,c) 



a=F u ■ a'| < 2. (6) 
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Here the equations £*(a) = u ■ a and C^(a') = u ■ a' 
indicated by Malus's law are used. To derive three-qubit 
Leggett's inequality, let's consider three 4-tuple settings 
(aj,a£, bj, Cj)(i = 1,2,3). The angle between all pairs 
(aj,a£) is the same 9 and (a' ; — a,) = 2sin|ei, where 
{ei,e2,e3} form an orthogonal basis. Then combining 
the three corresponding inequalities (^|) and the fact that 

Z)i=i l e « • u l > one obtains Yh=i |£* BC ( a i> b i) c i) + 
£* BC (a-, hi, Cj)| + 2 1 sin || < 6. Doing an integration over 
D(u, v, s) and using the fact that J dudvdsD(u, v, s) = 1, 
we obtain a three-qubit Leggett-type inequality: 




E 



Q«i + Qi 



<6, 



(7) 



where (2^ = J dKD(K)Cf c (a. ] ,h i ,c i )(aL j 



if J 



; and a. 



a if j 



i') is the physically measurable 
correlation function. Inequality (m) is our final three- 
qubit Leggett-type inequality. This inequality is a gener- 
alization of the two-quibt inequality obtained in Rcf . [[l3| . 
To obtain this inequality, we used the eight inequali- 
ties (|4a|-[lb|). In fact, there is an alternative approach, 
which is a little bit easier to derive three-qubit Leggett- 
type inequality, if we note the fact that the following 
equation is valid for any a, /3, 7 = ±1: 



I a ± Py\ =F a/37 = 1 



(8) 



Doing an integration J d\p u ^ VtS (X) for both side of 
Eq. (||), one can easily get inequality (j^) and hence in- 
equality (Q) by the same method. It is worthwhile to 
underline that the second approach does not require the 
eight inequalities coming from the positivity of probabil- 
ity distribution. It only relies on a very simple algebraic 
Eq. (g). Thus, it is, to some extent, a little bit easier. 

Quantum prediction for the correlation functions 
is Qijk = (a* • <7 ® bj ■ & (8 Cfc • a), where a 
arc Pauli matrices. For three qubits, we consider 
two types of entangled states: GHZ state and W 
state. We first look at the GHZ state \^)f HZ = 
l/\/2(|000) + I HI))- By choosing the measurement set- 
tings as ai = (cos |, — sin |, 0), a 2 = (0, cos |, — sin |), 



( 



a 3 

(0, cos |, sin |), a 3 



,0), 



,0), 



(sin |, cos ■ 



(cos f , 

2,^2/' "3 ~ v"" 2 ,v~« 2 ,0), bi = ci = (1,0,0), 

and b 2 = b,3 = C2 = C3 = -^,0), it is found that 
the quantum value of the left-hand side of inequality (Q) 
i sin I ) , which is larger than 6 when 
Moreover, noting that 6(cos | + 
= 2\/10sin(| + #0)7 the maximal quantum viola- 



Q _ 



< 



6(co S 2 , 3 „ 
4 arctan | 



is 

< e 

isinf 

tion is 2-\/l0 and it occurs when = 2 arctan | « 36.9°. 
The maximal violation of inequality (]7|) by the GHZ 
state is the same as that of the inequality given in 
Ref. |l3| by singlet state. We then consider the quan- 
tum violation of inequality (0) by generalized W states 
|^)W = sin £ cos 7? 1 100) + sin£sin?7|010) + cos£|001). In 
Fig. [TJ, we show the numerical results of quantum viola- 
tion by the family of generalized W states for the cases 



FIG. 1: (Color online) Numerical results of quantum viola- 
tions of inequality (Q) by the generalized W states = 
sin £ cos 77 1 100) + sin £ sin 77 1 010) +cos £|001). Quantum predic- 
tions are plotted versus the variation of 77 with £ = 7r/12, £ = 
tt/6, £ = 7r/4, £ = tt/3,£ = 5tt/12,£ = tt/2, and 77 is from to 
7r/2. Inequality (Q) is violated for the states whose quantum 
predictions are larger than 6. 



i = it/12, i = tt/6, £ = tt/4, £ = tt/3, £ = 5tt/12, £ = tt/2. 
It is clear that for some values of £ and 77, the inequality 
(0) is violated by l^)^. Furthermore, as pointed out in 
Ref. (TJJ , an arbitrary three-qubit pure state can be writ- 
ten in terms of five parameters as (up to local unitary 
transformations) : 



MulOOO) + 



/Mie^|100) 



/ti|101) 



(9) 



with /ij > 0, J^i^i = 1j an d < < 7r. For 1^)3, our 
numerical results show that the inequality (0) is also vio- 
lated for some region of the parameters and the maximal 
violation happens at = /J4 = | and /ii = /J2 = [i<x = 0, 
corresponding to the GHZ state. 

Our three-qubit Leggett-type inequality (Q) can be eas- 
ily generalized to the general n-qubit case. For all these 
inequalities, Alice has six measurement settings while 
other parties each has three settings. In general, the 
n-qubit Leggett-type inequality reads: 



E 



Qi 



<6, 



(10) 



where Qji...i arc physically measurable n-qubit correla- 
tion functions. The inequality ( fio| ) is violated for var- 
ious quantum entangled states. Its maximal violation 
happens for n-qubit GHZ state |"0)^ HZ = -^=([0 • • • 0) + 

1 1 - - - 1) ) and the maximal violation is 2vT0, too. In 
fact, for four-qubit case, we numerically calculated the 
quantum violations for all the pure entangled states and 
find that the maximal quantum violation for four-qubit 
Leggett-type inequality (|Io|) is indeed 2vTTj- 

To summarize, we have derived multipartite Leggett- 
type inequalities by using two different approaches. 
These inequalities are generalizations of the two-qubit 
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Leggett-type inequality obtained in Rcf. The max- 
imal violation of these inequalities is always 2^/10, and 
thus it does not scale with the system size. This shows 
a shortcome of our generalized inequalities. Of course, 
to acquire larger quantum violation and hence smaller 
experimental visibility threshold, one can think about 
increasing the number of measurement settings, as dis- 
cussed in Ref. jL3). However, increasing the number of 
settings will make the experimental tests more difficult. 
To obtain stronger multipartite Leggett-type inequalities 
with few settings is a very challenging problem. We 
hope that the two approaches presented in this paper 
will stimulate other works in this direction. Moreover, 
so far all the experimental tests of Leggett-type inequal- 
ities suffer from the "detection inefficient" loophole. A 
promising loophole free experimental scheme is to test 



these Leggett-type inequalities in topologically ordered 
systems that host non-Abelian anyons. Candidate sys- 
tems including quantum Hall liquid []l7|, rotating Bose 
condensates [Q, quantum spin systems [[Hi as well as 
Aharonov-Casher phase p3|. The entangled states can 
be obtained by braiding non-Abelian anyons ^J, while 
the measurements is implemented by fusing anyons to- 
gether. We will investigate this scheme subsequently. 

J.L.C. is supported by National Basic Research 
Program (973 Program) of China under Grant No. 
2012CB921900 and NSF of China (Grant Nos. 10975075 
and 11175089). This work is also partly supported by the 
IARPA MUSIQC program, the ARO and the AFOSR 
MURI program, and partly by National Research Foun- 
dation and Ministry of Education, Singapore (Grant No. 
WBS: R-710-000-008-271). 



[1] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 
777 (1935). 

[2] J. S. Bell, Physics (Long Island City, N.Y.) 1, 195 (1964). 
[3] J. Clauser, M. Home, A. Shimony, R. Holt, Phys. Rev. 

Lett. 23, 880 (1969). 
[4] D. Collins, N. Gisin, N. Linden, S. Massar, and S. 

Popescu, Phys. Rev. Lett. 88, 040404 (2002). 
[5] E. G. Cavalcanti, C. J. Foster, M.D. Reid, and P. D. 

Drummond, Phys. Rev. Lett. 99, 210405 (2007). 
[6] N. D. Mermin, Phys. Rev. Lett. 65, 1838 (1990); M. 

Ardehali, Phys. Rev. A 46, 5375 (1992); A. V. Belinskii 

and D. N. Klyshko, Phys. Usp. 36, 653 (1993). 
[7] J. L. Chen, C. F. Wu, L. C. Kwek, and C. H. Oh, Phys. 

Rev. Lett. 93, 140407 (2004). 
[8] S. J. Freedman and J. F. Clauser, Phys. Rev. Lett. 28, 

938 (1972); A. Aspect, P. Grangier, and G. Roger, Phys. 

Rev. Lett. 49, 91 (1982); Z. Y. Ou and L. Mandel, Phys. 

Rev. Lett. 61, 50 (1988); P. G. Kwiat, K. Mattle, H. We- 

infurter, A. Zeilinger, A. V. Sergienko and Y. Shih, Phys. 

Rev. Lett. 75, 4337 (1995); M. A. Rowe, D. Kielpinski 

, V. Meyer , C. A. Sackett , and W. M. Itano, Nature 

(London) 409, 791 (2001). 
[9] A. Cabello and J. A. Larsson, Phys. Rev. Lett. 98, 

220402 (2007); P. M. Pearle, Phys. Rev. D 2, 1418 (1970); 

E. Santos, Phys. Rev. A 46, 3646 (1992). 
[10] A. J. Leggett, Found. Phys. 33, 1469 (2003). 
[11] S. Groblacher, T. Paterek, R. Kaltenbaek, C. Brukner, 

M. Zukowski, M. Aspelmeyer, and A. Zeilinger, Nature 

(London) 446, 871 (2007). 
[12] T. Paterek, A. Fedrizzi, S. Groblacher, T. Jennewein, M. 



Zukowski, M. Aspelmeyer, and A. Zeilinger, Phys. Rev. 
Lett. 99, 210406 (2007). 
[13] C. Branciard, N. Brunner, N. Gisin, C. Kurtsiefer, A. 
Lamas-Linares, A. Ling, and V. Scarani, Nature Phys. 4, 
681 (2008). 

[14] C. Branciard, A. Ling, N. Gisin, C. Kurtsiefer, A. Lamas- 
Linares, and V. Scarani, Phys. Rev. Lett. 99, 210407 
(2007) 

[15] W. Laskowski, T. Paterek, M. Zukowski, and C. Brukner, 

Phys. Rev. Lett. 93, 200401 (2004). 
[16] A. Acm, A. Andrianov, L. Costa, E. Jane, J. I. Latorre, 

and R. Tarrach, Phys. Rev. Lett. 85, 1560 (2000); A. 

Acm, A. Andrianov, E. Jane, and R. Tarrach, J. Phys. 

A 34, 6725 (2001); C. Emary and C.W. J. Beenakker, 

Phys. Rev. A 69, 032317 (2004). 
[17] J. S. Xia, W. Pan, C. L. Vicente, E. D. Adams, N. S. 

Sullivan, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer, K. W. 

Baldwin, and K. W. West, Phys. Rev. Lett. 93, 176809 

(2004). 

[18] N. R. Cooper, N.K. Wilkin, and J. M. F. Gunn, Phys. 

Rev. Lett. 87, 120405 (2001). 
[19] M. Freedman, C. Nayak, K. Shtengel, K. Walker, Z. H. 

Wang, Ann. Phys. (N.Y.) 310, 428 (2004); P. Fendley 

and E. Fradkin, Phys. Rev. B 72, 024412 (2005). 
[20] J. L. Chen, H. Y. Su, C. F. Wu, D. L. Deng, and C. H. 

Oh, arXiv: 1104.4173 (2011). 
[21] D. L. Deng, C. F. Wu, J. L. Chen,and C. H. Oh, Phys. 

Rev. Lett. 105, 060402 (2010). 



